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One of the schemes that has been proposed for the adaptive control of large 
flexible space structures is shown in figure 1. This approach is based upon a 
modal decomposition of the dynamic response of the flexible structure and is 
designed to make use of the parallel processing features of modern minicomputers. 
Satisfactory performance of the parallel structure identification technique shown 
in figure 1 can he achieved only when the approximation functions noted in the 
figure correspond to the natural modes of the flexible structure. The work 
su mma rized here presents a technique for estimating both mode shapes and modal 
parameters. 
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(a) Motion of flexible structure. 


q ( k + I ) = A, q ( k ) + q ( k - I ) + B } U ( k ) + B 2 U ( k - I ) 


y ( k ) = H ( k ) + w ( k ) 

q 

(b) Modal description. 


Figure 1.- Problem formulation. 
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Figure 2 shows the analytic background for the formulation of the on-line 
identification problem. The motion of the flexible structure is expressed as a 
sum of NM terms involving the modal shape function £^(s) and the modal amplitudes 
q^Ctjp^I where s denotes the spatial variable and t represents time. The 
parameter vector p£ represents a set of parameters l^l±y^2±9^1±9^2±^ f° r ea ^h 
mode. The identification problem is formulated in terms of the modal description 
shown in figure 2 (B), where the Ai,i=l ,2 are diagonal matrices of order NM, U(k) 
denotes a vector of actuator signals and the B-^ are rectangular matrices. Matrix 
H, relating the modal amplitudes to the measurement vector y(k), has columns that 
are linearly independent in the modal model. On-line measurements {y(k) } , {U(k) } 
are to be processed in order to obtain estimates of the matrices A^, B^, and H. 




Figure 2 .- Distributed adaptive control. 
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Algorithms for on-line parameter identification are shown in figure 3. 
Assuming that the approximation functions are given, the modal parameters are 
updated using an output error formulation where the error e^ is given by 


e^k-1) = qi (k-l) - A liqi (k-2) - A 2±q± (k-3) - S li F i (k-2) - B^Ck-S) 


where Fi(k) 
satisfy 


denotes a modal force component. 


The weights 


W. 


are selected to 


W l q i (k-2) + W 2 q i (k-3 > + W 3 F 2 (k-2) + W 4 F 2 (k-3) < 2 


to insure stability of the identification algorithm. 

The approach for updating the approximation functions is based on combining 
the modal equations shown in figure 2 into the single relation involving matrices 
M^ and where 


M ± = HA ± H 


N ± = HB ± 


i = 1,2 


and (H) + denotes the pseudo-inverse of H. Using a regression analysis approach, 
data are collected over a time interval of length N and a least-squares estimate 
of M-j^ and Ni is obtained. When the number of modes in the modal approximation 
is the same as the number of sensors, the eigenvectors of Mi correspond to the 
columns of H. 
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(a) Pole-zero characteristics - output error formulation. 


y ( k + I ) = M | y(k)+M 2 y(k-l) + N ( U(k) + N 2 U(k-ll + n(k + l) 


N £ Y < N ) = S(NIM T + V(N); m t = 


« T = S + ( N > Y ( N ) 


(b) Approximation functions - regression analysis. 


Figure 3.- On-line parameter identification. 
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A schematic diagram of the flexible beam used to test the identification 
algorithms is shown in figure 4. Sensor and actuator locations along the beam are 
shown along with the modal frequencies and mode shapes for the first three flexible 
modes obtained from the SPAR analysis program. 
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Figure 4.- Modes of interest. 







Figure 5 shows the modal decomposition for mode 3 unforced response. A set 
of approx ima tion functions, obtained from the SPAR analysis was used and eight 
modes were assumed in the real-time program that produced the modal amplitude 
signals shown. For this response there is some excitation of modes 1 and 2. 
However, higher frequency modes do not appear to be excited. 
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Figure 5.- Modal decomposition. 
Eight SPAR modes; mode 3 excited. 
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Figure 6.- Output of par 







figure 7 shows the outputs of four on-line identifiers resulting from a 5-Hz 
excitation produced by actuator 1. Note that the parameters of the first vibration 
mode are identified b.ut the high, frequency modal parameters are not found as a 
result of the 5-Hz excitation. 
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Figure 7.- Output of parameter identifiers. SPAR derived values 
Sinusoidal excitation by actuator 1 at 5 Hz . 
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f Figure 8 demonstrates the result of applying the mode shape identification 

procedure (a) to a simulation wherein the beam was given an arbitrary displacement 
and (b) in an experiment wherein only the first flexible mode was excited. In the 
simulation the normalized mode shape obtained agreed with that produced by the 
SPAR analysis. However, in the laboratory experiment a distortion in the identified 
mode shape was observed. In both cases the mode shape was identified on line after 
approximately 1 second of data processing. 



(b) Experimental 
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Figure 8.- Identification of mode shape. 
First flexible mode. 
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A table of identified modal parameter values that resulted from a simulation 
incorporating both the mode shape and modal parameter identification schemes is 
shown in figure 9. Identification of mode shapes and parameters was accomplished 
in approximately 1 second of processing measurements from eight sensors distributed 
along the beam* 


[Simulation - all modes excited at parameter values] 


Flexible mode 

Initial 

Final 

Actual 

i 

-0.8495 

1.8726 

1.87267 

2 

-1.7194 

1.0991 

1.099181 

3 

0.09918 

-0.71936 

-0.71940 

4 

0.83907 

-1.9960 

-1.9960 

5 

0.8726 

0.15047 

0.15049 

6 

-2.996 

1.8387 

1.83907 


Figure 9,- Identification of modal parameters* 


SUMMARY 


• Presented algorithms for on-line parameter and mode-shape estimation 

• Examined identification performance using computer simulations and a 

limited number of laboratory experiments 

• Future work will include: 

(1) Further experimental studies 

(2) Development of separable nonlinear least-squares approach 

to identification 

(3) Development of on-line performance measures for sequential 

processing decisions 
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